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ABSTRACT 



We find that the recently reported departure from the Eddington luminosity limit for the 
highest quasar black hole masses at a given redshift is an artifact due to biases in black hole 
mass measurements. This sub-Eddington boundary (with non-unity slope) in the quasar mass- 
luminosity plane was initially reported by Steinhardt & Elvis (2010a) using the FWHM-based 
black hole mass catalogue of Shen et al. (2008). However, the significance of the boundary is 
reduced when the FWHM-based mass-scaling relationship is recalibrated following Wang et 
al. (2009) and using the most updated reverberation mapping estimates of black hole masses. 
Furthermore, this boundary is not seen using mass estimates based on the line dispersion of 
the same quasars' Mgll emission lines. Thus, the initial report of a sub-Eddington boundary 
with non-unity slope was due to biases in estimating masses using the FWHM of a fit of 
one or two Gaussians to quasar Mg II emission lines. We provide evidence that using the line 
dispersion of the Mg II line produces less biased black hole mass estimates. 
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1 INTRODUCTION 

Steinhardt & Elvis (2010a) (hereafter SElOa; see also Steinhardt 
& Elvis 2010bcd) have claimed that there exists a departure from 
the Eddington luminosity (LEdd) boundary for the highest quasar 
black hole masses at a given redshift; in other words, a sub- 
Eddington boundary (SEB) with non-unity slope. SElOa investi- 
gated the authenticity of this SEB and argued that its presence was 
neither due to measurement error nor due to small number statis- 
tics. 

SElOa used the black hole (BH) masses of Shen et al. (2008), 
who adopted mass scaling relations for H/3 and C IV from Vester- 
gaard & Peterson (2006) and for Mg II from McLure & Dunlop 
(2004). Both those references have made two key assumptions. 
They have assumed a tight relationship between the radius of the 
emitting region and the AGN continuum luminosity in the form of 
R oc L'^, where /3 ~ 0.5. They have also assumed a relationship 
between the characteristic virial velocity of the gas emitting a par- 
ticular emission line, Vviriai, and the single-epoch FWHM of the 
same line of the form v^iriai oc FWHM. SElOa argued that depar- 
tures from the above assumptions sufficient to affect the presence of 
a SEB were unlikely, although they acknowledged that the slope of 
a SEB is dependent upon the BH mass scaling relation used (their 
§4.5.2). 

Peterson et al. (2004) suggested that the BLR virial veloc- 
ity can be characterized by either the line dispersion (JKne or the 
FWHM measured directly from the emission line profile in an rms 



spectrurrQ obtained during a reverberation mapping study (see, 
e.g., Figure 3 of Peterson et al. 2004). This result has been inter- 
preted as arising from linear relations between virial velocity and 
each of the rms FWHM and rms aune- The same study shows 
that the line dispersion measured from a mean spectrum is also 
linearly proportional to the virial velocity; however, the FWHM 
measured from a mean spectrum may not be linearly proportional 
to the virial velocity (see e.g.. Figure 4 & 5 of Peterson et al. 
2004). Since single-epoch spectra are more similar to a mean spec- 
trum than an rms spectrum, Wang et al. (2009) claim that the re- 
lation Vviriai OC FWHM may not hold in general when we use 
single epoch spectra and suggest using a relationship of the form 
""virial OC FWHM^ instead. The virial relation still holds the same 
form of XbH °^ ^ '"liriai- 

Motivated by Wang et al. (2009), we investigate the reliabil- 
ity and biases of the SEB by re-estimating the mass scaling rela- 
tion and by using two different methods of estimating line width; 
namely, using the FWHM and the line dispersion of a single epoch 
spectrum. We present and discuss the SEB in §(2] We compare the 
use of FWHM- and line dispersion-based BH masses in §[3] We 



In a reverberation mapping study, the mean spectrum is defined as 



F(X) = X/i-i where Fi{\) is the (th spectrum of the N spec- 

tra that compose the reverberation data for one object. The rms spectrum is 

defined as 5(A) = { X^ili [■^»(-'^) ~ ^WfY^^ (Peterson et al. 
2004). 
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discuss the effects of using different mass scaling relationships on 
the SEB in §|4] We end with our conclusions in §|5] 



2 THE SUB-EDDINGTON BOUNDARY 

SElOa have plotted the mass-luminosity plane using the black hole 
masses of Shen et al. (2008) from the SDSS DR5 quasar sample 
(Schneider et al. 2007). They have used 62185 quasars over the 
redshift range 0.2 < z < 4.0 with three BH mass scaling relations 
using the FWHMs of the H/3, Mg II, and C IV emission lines. 

SElOa divide their H/3, Mgll, and CiV mass estimates into 
13 redshift bins (though their conclusions are only based on H/3 
and Mg II). They show that SEB of non-unity slope is a prominent 
feature in most of the bins, especially the highest redshift ones. The 
Mg II emission line can be measured from SDSS spectra in eight of 
these redshift bins. 

In Figure [T^, we illustrate the SEB seen by SElOa using the 
Shen et al. (2008) quasar BH masses in a particular redshift bin. The 
gap between the red dot-dashed line (showing an Eddington ratio 
of one) and the black dot-dashed curve (tracing the 95th percentile 
of the observed distribution) is called the sub-Eddington boundary. 

In Figure |2] panels (a) to (h), we use the same quasar sample 
as Figure [T^ but plot the mass-luminosity plane for 8 redshift bins 
between 0.76 < z < 1.98. Based on a similar diagram (their Fig- 
ure 8, panels 4 to 11), SElOa concluded that the most luminous 
low-mass quasars in any given redshift bin are at the Eddington 
limit but that the most luminous high-mass quasars are not at the 
Eddington limit. 

However, a SEB with non-unity slope cannot be seen using the 
Rafiee & Hall (201 1) BH mass catalogue which is based on the line 
dispersion of Mg II instead of its FWHM. This catalogue consists 
of 27602 quasars with 0.7 < z < 2.0, where the redshift limits are 
set by our method of removing the Fe II pseudocontinuum from the 
Mgll region of the spectrum. In Figure [TJ) and Figure [3] panels (a) 
to (h), we plot the mass-luminosity plane for the same objects as 
in Figure|2]but using the Rafiee & Hall (201 1) BH mass estimates. 
This plot shows no evidence for the existence of an SEB with slope 
different from unity. Instead, in the lower redshift bins there is a 
sub-Eddington boundary whose slope is consistent with unity; i.e., 
a boundary at a fixed fraction of the Eddington limit. As the redshift 
of the bin increases, the boundary moves to a higher fraction of the 
Eddington limit until it reaches the Eddington limit itself. 

Comparing Figure [T^ with Figure [T}3 or comparing Figure |2] 
with Figure [3] shows that two BH mass samples from the same 
SDSS data set yield different distributions in the mass-luminosity 
plane. The fact that different mass-scaling relations were used for 
these samples suggests that the SEB is sensitive to the mass cali- 
bration or the line width indicator used, or both. We now proceed to 
investigate these line width indicator and mass calibration issues. 



3 CHARACTERIZING THE LINE WIDTH 

Shen et al. (2008) fit Mg II as a sum of two Gaussians, one con- 
strained to have FWHM<1200 km s"\ and use the FWHM of the 
broader Gaussian as the Mg II FWHM. It is convenient to assume 
that quasar emission lines can be modeled as sums of two Gaus- 
sians; however, the validity of this Gaussianity assumption has not 
yet been fully investigated. There are only a few studies on the 
non-Gaussianity of quasar broad emission lines in the context of 
BH mass estimation. Peterson et al. (2004) have suggested using 



the line dispersion and not the FWHM of the emission profile for 
mass measurements for several reasons including the Gaussianity 
assumption. Collin et al. (2006) have used both the FWHM and 
rms line dispersion of reverberation-mapped AGN to investigate 
the Gaussianity of the lines. They reported two classes of objects 
according to the value of FWHM/crHne: Pop. I with values be- 
low the Gaussian ratio of 2\/2 In 2 ~ 2.35 and Pop. II with values 
above it. Peterson (2007) have also reported that for H/3 this ratio 
spans 0.71 < FWHM/aane < 3.45 with an average ratio of 2.0, 
indicating that the emission lines may not be well fit by single or 
double Gaussians. 

We have investigated the significance of the Gaussianity as- 
sumption for BH mass estimates using the line dispersion aune of 
all objects from Rafiee & Hall (201 1) and the FWHM reported by 
Shen et al. (2008) for the same objects. In Figure|4l we have plot- 
ted the histogram of the resulting Mg II FWHM / aune distribution. 
We report a range of 1 < FWHM/ aune < 5 with a mean value 
of 2.7 and a mode of 2.55 (the blue dash-dotted line in Figure|4]l. 

In Figures |5^ and [Sj), we show contour plots of the distri- 
bution of these two Mgll line width indicators as a function of 
their ratio. There is a strong correlation of FWHM with the ra- 
tio FWHM/(j;i,ie, but no significant correlation for anne- Figure 
[5}; plots FWHM versus o";i„e- This contour plot is another way of 
demonstrating that the distribution of FWHMJauns (e.g., between 
the five reference lines plotted) as a function of aune shows only 
moderate variation from an average distribution, while as a function 
of FWHM the distribution of FWHM/crHne shows great variation. 
Recall that Peterson et al. (2004) found that while the virial velocity 
in the BH mass relation is well characterized by aun^ , it is less well 
characterized by FWHM (§1). Combining that result and ours, we 
can say that aii„e is a good representation of the virial velocity re- 
gardless of the line shape (as crudely measured by FHWM/aune), 
whereas FWHM is not. This change in the typical ratio between 
FWHM and anne with FWHM again calls into question the use of 
a BH mass relation calibrated to the FWHM. 

In Figure |6l we show contour plots of quasar BH masses ver- 
sus FWHM/(THne for three scenarios. Panel (a) shows the Shen et 
al. (2008) results which assume a Gaussian profile for Mg II, and 
panel (b) the Rafiee & Hall (2011) results which directly use the 
line dispersion. Figure[6^ shows that the dependence of FWHM on 
the FWHM/(7ii„e ratio means that a broader range of BH masses 
are found using FWHM-based estimation as compared to aun^- 
based estimation (Figure |6j5). Figure |6j; shows that the recalibra- 
tion of FWHM-based masses discussed in the next section helps to 
reduce the discrepancy between the range of BH masses estimated 
using FWHMs and that estimated using aune values. 



4 RECALIBRATING THE BLACK HOLE MASS 
SCALING RELATION 

4.1 Statistical Methods 

We have used numerous different methods to estimate the best-fit 
parameters of two- and three-parameter BH mass scaling relation- 
ships. However, we base our conclusions on the MLINMIX_ERR 
method (Kelly et al. 2007), also used by Wang et al. (2009), which 
we report in Table[T] The MLINMIX_ERR method is preferred over 
other methods because it is the only fitting method which takes in- 
trinsic scatter and uncertainties in both parameters into account. It 
should be the default fitting method in such cases, though its results 
may not differ significantly from other methods for two-parameter 
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fits. We also use another Metropolis-Hastings Markov chain Monte 
Carlo simulation (MCMC; Haario et al. 2006) which gives very 
similar results to LINMIX_ERR and MLINMIX^RR; they both 
use the Bayesian approach to the regression problem but use differ- 
ent sampling algorithms. The MCMC and MLINMIX.ERR (and 
LINMIX_ERR) simulations yield confidence levels as well as po- 
tential outliers and the distribution of the acceptable parameters and 
their associated errors within the parameter space. 

Other methods used in previous studies are less sophisticated. 
For example, the BCES method (Akritas & Bershady 1996, e.g., 
used by McLure & Dunlop 2004, Kaspi et al. 2005, Vestergaard 
et al. 2006) considers bivariate correlated errors and a possible in- 
trinsic scatter but it does not report that scatter. A quantified value 
of the intrinsic scatter is available from the FITexy-T02 method 
(Tremaine et al. 2002, e.g., used by Kaspi et al. 2005, Vestergaard 
et al. 2006). However, FITexy-T02 does not account for the bivari- 
ate correlated errors in the fitting process. Nonetheless, for con- 
venience and to illustrate the large dispersion which results from 
applying non-optimal fitting algorithms, we report the results from 
all these methods in Table |2] 



4.2 The Wang et al. Three-Parameter Relationship 

Wang et al. (2009) have studied the mass scaling relations for the 
H/3 and Mg II emission lines in 29 out of 35 low redshift AGNs for 
which Peterson et al. (2004) have reported mass estimates made 
through reverberation mapping (RM) (see Table [TJ. Because the 
single-epoch FWHM of a given line is unlikely to be an exact tracer 
of its virial line width, in determining new mass-scaling relations 
they have assumed a more flexible relation between the virial line 
width and the FWHM in the form of vl^^^^i oc FWHM^ where 7 
is not fixed at 2, as in conventional virial relations, but rather is a 
free parameter. 

In their mass scaling, Wang et al. (2009) adopted the RM BH 
masses given in the Peterson et al. (2004) for 24 objects but used 
the RM BH masses from more recent reverberation campaigns for 
the other 5 objects: NGC 4593 from Denney et al. (2006), NGC 
4151 from Metzroth et al. (2006), PG 2130-1-099 from Grier et al. 
(2008), and NGC 4051 from Denney et al. (2009). They (and we) 
take the (j;i„e(H/3,rms) values of all 29 RM objects from Peterson 
et al. (2004). 

For their full 29-quasar RM sample, Wang et al. (2009) 
found a different power-law relation between FWHM(Mgll) and 

o-iine(H^,rms): 



log 



(7K„e(H/3,rms) 



1000 km s" 



= (0.85 ±0.21) log 



FWHM(Mg II) 
1000 km s"^ 
-(0.21 ±0.12). (1) 



Based on this result, Wang et al. (2009) concluded that the line- 
emitting locations of Mg II are different from those of H/3 in the 
broad line region (BLR). That conclusion may be premature, as the 
difference from unity slope is < la statistical significance and the 
FWHM(Mg II) values are measured from single epoch spectra and 
not from the same rms spectra as the ounc(iil3) values. Nonethe- 
less, they suggest fitting a three parameter relation between Algjj, 
FWHM and AL3000: 



log 



a^bh(-R^) 



a ± /3 X log 
+7 X log 



1044 erg s-i 
FWHM(Mg II) 



1000 km s" 



(2) 



where 7 = 1.7 ± 0.42 provides consistency between Mgll-based 
single-epoch BH masses and H/3-based RM BH masses. 



4.3 An Up-to-date, Recalibrated Three-Parameter 
Relationship 

Motivated by Wang et al. (2009), we have revisited the mass scaling 
relation used by Shen et al. (2008), which was derived by McLure 
& Dunlop (2004) based on 20 RM objects studied in McLure & 
Jarvis (2002). McLure & Dunlop (2004) reported a one-to-one 
relation between FWHM-based TWgjj (H/3) and AlgH 
logA1gH(H/3) = 1.00(±0.08) logXBH(Mgii)±0.06(±0.46) 
and concluded that H/3 can be replaced by Mg II for purposes of 
mass estimation. 

We have updated the BH mass list used by Wang et al. (2009) 
with the most recent RM BH mass estimates of NGC 5548 (Bentz 
et al. 2009) and NGC 3227 (Denney et al. 2010). Using those up- 
dated BH masses, we have recalibrated the Wang et al. (2009) mass 
relation. Comparing results before and after updating the RM sub- 
sample shows no statistically significant change in the fitting re- 
sults, except in those of FITexy and some improvement in the value 
of the correlation coefficient in the three-parameter fitting results. 
This improvement means that the regression line fits the updated 
sample better than that of Wang et al. (2009) or McLure & Jarvis 
(2002). 

Here, we report the three parameter mass scaling relations us- 
ing the most updated RM masses and the MLINMIX_ERR regres- 
sion package. 



log 



Mbh(-R^) 



IO6M0 



(1.25 ±0.22) ±(0.51 ±0.08) 



X log 



1044 erg s- 



±(1.27 ±0.40) X log 



±a 



/ FWHM(Mgll) 



1000 km s" 



±(0.15 ± 0.5) dex(intrinsic scatter). 



(3) 



where '^iag[Jv{-QY{/w^MQ\^^^^^^^^^'^^^ statistical error of the 
black hole mass from: 



[0.048 ± 0.0012 (In (£))^ ± 0.048 
±0.030 (In (o;))^ ±0.32^ 



ae 
P 



(4) 



where £ — [AL3000/IO*''] is in units of erg s ^, cj — 
[FWHM/1000] is in units of km s"\ and and a I are estimated 
errors from our fitting process. Like Wang et al. (2009), we find 
that the relationship of A^bh '° ^8 ^' FWHM follows a power law 
with a smaller value than 2 but that the deviation is not statistically 
significant (< 2(t). 



4.4 Recalibrated Mass-Luminosity Plane 

We have shown in Figure [3] that when the Rafiee & Hall (2011) 
sample shows a SEB in the mass-luminosity plane, its slope paral- 
lels that of the Eddington limit, unlike what was found by SElOa. 
The same is true with most of the new mass scaling relations. The 
mass-luminosity plane using the Shen et al. (2008) measurements 
with a re-scaled mass-relation assuming oc FWHM^ with 



4 A. Rafiee and P. B. Hall 



7 < 2 shows no sign of a SEB with non-unity slope (see FigureQ. 
However, the slope of the SEB may be less than unity if 7 ^ 2; as 
a change from 7 = 2 to 7 < 2 has increased the slope of the SEB, 
a change from 7 = 2 to 7 > 2 should decrease its slope. 

Our results show that a recalibration of the mass scaling re- 
lation can produce different shifts in the high and low mass tails 
of the distribution, rather than a uniform overall shift to lower or 
higher masses. Compared to masses recalibrated with the latest 
mass-scaling relation, the published Shen et al. (2008) masses over- 
estimate the highest BH masses and underestimate the lowest BH 
masses. The overall distribution of the log difference of mass is 
nearly symmetric around zero, as shown in Figure [8^. However, 
the new mass estimates versus the Shen et al. (2008) masses de- 
picted in Figure [8j) show a rotation in the mass distribution around 
log A^gjj ~ 9 instead of a systematic shift in the new mass es- 
timates. This rotation around log A^gj^ ~ 9 (Figure[8j>) may ex- 
plain the presence of a SEB of non-unity slope in the original Shen 
et al. (2008) data, especially at the highest masses. 

Another potential cause of differences between BH mass 
datasets is the value of P used for the luminosity dependence in 
the mass calibration equation. The SElOa mass scaling relations in 
H/3, Mgll and CiV used /3 = 0.61, 0.62 and 0.53, respectively, 
whereas for Mgll we use /3 = 0.5 in Rafiee & Hall (2011) and 
/3 = 0.51 for the recalibrated three-parameter relationship above. 
Even assuming a perfect correlation between FWHM and aane, 
for two quasars of luminosity differing by a factor of 10 for which 
our method measures the same mass, SElOa would obtain a mass 
higher by 0.11 or 0.12 in the log for the higher luminosity quasar. 
That systematic difference will decrease the slope of the upper en- 
velope of the quasar distribution in the mass-luminosity plane, al- 
though not by enough to explain the full deviation from unity of the 
SEB slope in SElOa. 

4.5 Other Surveys 

Other surveys which use the FWHM as a line width indicator may 
yield biased BH masses if 7 = 2 is assumed. 

Kollmeier et al. (2006) use the FWHMs of H/3, Mg II and C IV 
measured directly from smoothed spectra from the AGES survey 
for their BH mass measurements at 0.3 < z < 4, with 7 = 2 and 
/3 = 0.88 for Mgll at 0.4 < 2 < 2. They do not report a SEB 
of non-unity slope, but may see a weak one in their 1 < z < 2 
redshift bin (their Figure 10, lowest three panels). 

Netzer & Trakhtenbrot (2007) use the FWHM of the sum of 
two broad Gaussians fit (along with one narrow Gaussian) to H/3 in 
SDSS spectra for their BH mass measurements al z < 0.75. They 
use a mass scaling relation with 7 = 2 but /3 — 0.65, and appear 
to find a SEB of non-unity slope (their Figure 1, bottom panel, and 
Figure 3, red and blue histograms). 

The high value of /3 in the above calibrations may partially 
compensate for any biases due to the use of FWHM or its calibra- 
tion. The monochromatic luminosity is little affected by the shape 
of the line (see Figure[5jl), while the FWHM clearly depends on the 
shape (see e.g., Figure|5^). A heavy weight given to the luminosity 
in the scaling relationship reduces the effect of the sensitivity of 
the FWHM to the line shape. This may reduces the presence of the 
SEB in, e.g., Kollmeier et al. (2006). 

4.6 Results using the SDSS Seventh Data Release 

Shen et al. (2010) have reported estimates for BH masses of 105783 
quasars in the SDSS Seventh Data Release (DR7). Their fiducial 



("SIO") Mgll-based BH mass estimates use the FWHM of the 
broad component of Mg II as fitted by up to 3 Gaussians in con- 
junction with a narrow Gaussian of width ^1200 km s~^. These 
mass estimates are calibrated by adopting the McLure & Dunlop 
(2004) slope /3 — 0.62 and setting the normalization a to match 
the H/3-based masses of Vestergaard & Peterson (2006) on aver- 
age. Figure [8}; shows that the BH mass estimates from Shen et 
al. (2010) have not significantly changed with respect to the mass 
catalogue reported by Shen et al. (2008). In Figure |9l we show 
the results of using this larger sample with updated Mgll-based 
mass estimates. Furthermore, following Wang et al. (2009), we can 
use the recalibrated mass-scaling relation reported in Equation [3] 
to recalculate the BH masses of objects within the redshift range 
0.76 < z < 1.98. The results, shown in Figure [Tol show no SEB 
of any kind in the mass-luminosity plane. Apparently, the presence 
of more low luminosity quasars in DR7 along with the correction 
applied to the mass-scaling relation have removed the signature of 
the SEB from the mass-luminosity plane. 



5 CONCLUSION 

SElOa have argued that the presence of a sub-Eddington bound- 
ary (SEB) with non-unity slope in the mass-luminosity plane con- 
structed from the Shen et al. (2008) FWHM-based BH mass data 
set is not due to measurement error and is statistically significant. 
They have shown that the SEB they found is not an artifact due to 
small number statistics. We have investigated the existence of such 
a SEB using our own data set (Rafiee & Hall 2011) which uses 
the line dispersion of Mg II to estimate BH masses. Furthermore, 
motivated by Wang et al. (2009), we have considered the effects of 
recalibrating the mass scaling relationship using the general three- 
parameter relation proposed in Equation [2] and taking into account 
the latest mass updates from reverberation campaigns. Finally, we 
have investigated the SEB using the larger SDSS DR7 BH mass 
catalogue of Shen et al. (2010). Overall, we conclude that: 

• We have found no sign of a SEB with non-unity slope in the 
Rafiee & Hall (2011) data set, before or after recalibration, nor in 
the recalibrated Shen et al. (2008) or Shen et al. (2010) data sets. 

• The non-unity slope of the SEB in the original data of Shen 
et al. (2008) arose from the mass-scaling relation used and is likely 
not a real feature in the quasar mass-luminosity plane. However, the 
implications of a sub-Eddington boundary of unity slope remains a 
worthwhile topic for theoretical investigation. 

• The presence or absence of a SEB of non-unity slope is rooted 
in the choice of line width indicator and in the calibration of the 
mass scaling relation used. Further studies of single epoch spectra 
of RM sample AGN may give clues about the roots of the problem 
with different line width indicators and calibrations. Until then, we 
suggest using the line dispersion, which is a more robust parameter 
to represent the line width, to estimate the BH masses of quasars 
when using a scaling relation of the form v^iriai oc aune- The 
FWHM can be used; however, the scaling relation may have a dif- 
ferent form vli^i^i oc FWHM^, with 7 < 2. 
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Figure 1. The quasar Mg II mass-luminosity plane, (a) using the original BH masses of Shen et al. (2008), and (b) using the BH masses of Rafiee & Hall 
(201 1). The lower horizontal dot-dashed line is the approximate lower luminosity limit for a low-redshift object at the faint magnitude limit of the SDSS. The 
upper horizontal dot-dashed line is the upper luminosity limit for a high-redshift object at the bright saturation limit of the SDSS. Objects below the detection 
limit or above the saturation limit are not in the SDSS spectroscopic sample. The vertical dot-dashed line represents an upper mass limit of logA^gjj = 10. 
SElOa have claimed there exists a sub-Eddington boundary of non-unity slope, which is here illustrated in panel (a) for objects with logA^gj^ > 8.2 and 
log L(,oi > 46.1 as the region between the black dot-dashed curve (the 95% upper luminosity limit of the distribution as a function of mass) and the red 
dot-dashed line (an Eddington ratio of one). 
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Figure 2. The quasar- Mg II mass-luminosity plane for 8 redshift bins using the original BH masses from Shen et al. (2008). A sub-Eddington boundary of 
non-unity slope for objects with the highest A^gj^ in a given redshift bin can be seen, as found by SElOa. 
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Figure 6. Contour plot of the distribution of log Mg II black hole mass versus FWHM/(T;.i„j, ratio, a) Using Shen et al. (2008) Mg II BH mass estimates (before 
recalibration). b) Using Rafiee & Hall (201 1) Mg II mass estimates, c) Using the FWHM from Shen et al. (2008) to estimate the BH masses via Equation[3] In 
all plots, the FWHM/au^^ ratio uses FWHM from Shen et al. (2008) and an^^ from Rafiee & Hall (2011). 



Biases in the Quasar Mass-Luminosity Plane 9 



7 48 

CO 

D) 47 



0.76<z<0.912 ^ 



CD 



C3) 
O 



46 
45 




7 48 
CO 

O) 47 



0.912 <z< 1.06 



CD 



O 



46 
45 



10 



LogMBH^gii(Shen08) 



1.22<z< 1.37 




CD 



o 



LogMBH^gii(Shen08) 



1.68 < z < 1.83 




7 48 

CO 

D) 47 



1.06 <z< 1.22 



CD 



O 




10 



10 



LogMBH^gii(Shen08) 



LogMBH^gii(Shen08) 



7 48 
CO 

O) 47 



1.37 <z< 1.52 



46 
45 



(e) 



7 48 

C/3 

O) 47 



1.52 <z< 1.68 



CD 



O 




10 



10 



LogMBH^g,i(Shen08) 



LogMBH^gii(Shen08) 



1.83 <z < 1.98 




LogMBH^gii(Shen08) 



LogMBH^gii(Shen08) 



Figure 7. The quasar Mg II mass-luminosity plane for 8 redshift bins when we apply a new mass-scaling relation to the original FWHM and luminosity of 
Shen et al. (2008). The slope of the upper envelope of the quasar distribution is consistent with unity in all redshift bins. 



Table 2. Re-scaling results using most updated black hole masses. 
logCA^BH/lO"-^©) - CSlogCALA/lO^^ergs-i) = a -I- 7 X 
log(FVKf/A//1000kms-i). 



Fitting Metliod 


Intercept(a) 


Slope(7) 


Intrinsic Scatter 


OLS(Y|X) 

OLS(X Y) 

OLS Bisector 

OLS Bisector bootstrap 

OLS Bisector Jacknife 


1.244 ± 0.25 
-0.136 ±0.38 
0.85 ± 0.24 
0.833 ± 0.26 
0.849 ± 0.29 


1.27 ± 0.47 
3.92 ± 0.66 
2.03 ± 0.42 
2.05 ± 0.44 
2.03 ± 0.50 




FITexy 
FITexy-T02 


0.728 ± 0.09 
0.85 ± 0.23 


1.99 ±0.17 
1.95 ±0.43 



0.29 


BCES(Y|X) 

BCES(Y X) bootstrap 

BCES(X Y) 

BCES(X Y) bootstrap 

BCES Bisector 

BCES Bisector bootstrap 

BCES Orthogonal 

BCES Orthogonal bootstrap 


1.21 ± 0.27 
1.16 ±0.33 
0.0343 ± 0.34 
0.203 ± 202.00 
0.846 ± 0.25 
0.822 ± 0.28 
0.259 ± 0.31 
0.346 ± 100.00 


1.34 ±0.51 
1.41 ±0.60 
3.95 ± 0.58 

3.15 ± 475.00 
2.03 ± 0.44 
2.07 ± 0.49 
3.16 ±0.51 

2.93 ± 235.00 




MCMC (median posterior distribution) 
MCMC (mean posterior distribution) 


1.24 ±0.20 
L24±0.20 


1.28 ± 0.36 
1.27 ± 0.36 


0.37 ± 0.05 
0.37 ± 0.05 


LINMIX_ERR (median posterior distribution) 


1.25 ±0.22 


1.21 ± 0.40 


0.11 ±0.05 
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Figure 8. a) The log difference between the Shen et al. (2008) BH masses using the old and new calibrations. There is no systematic shift in the mass 
distribution, b) Scatter plot of the same BH masses. The mass distribution has rotated around the value log A^gfj ~ 9 under the new calibration, causing the 
SEB to disappear, c) Contour plot of BH mass estimates from Shen et al. (2010) versus Shen et al. (2008) for the Mg II redshift range. 
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Figure 9. The quasar Mg II mass-luminosity plane for 8 redshift bins using tlie fiducial "S 10" BH masses from Shen et al. (2010). A sub-Eddington boundary 
of non-unity slope for objects with the highest A^gjj in a given redshift bin can be seen partially in panels (a) to (e), as found by SElOa. 
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Figure 10. The quasar Mgll mass-luminosity plane for 8 redshift bins using Equation[3]to re-estimate the "SIO" BH masses from Shen et al. (2010). No 
sub-Eddington boundary of any form, with non-unity or unity slope, can be seen here. 



